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1. Introduction 

We suppose that a discrete eigenvalue E n — F n (v) of the Schrodinger 
Hamiltonian 

H = -A + vf(x) 

is known for all sufficiently large values of the coupling parameter v > v c 
and we use this data to reconstruct the potential shape f(x). The usual 
'forward' problem would be: given the potential (shape) f(x), find the cor- 
responding energy trajectories F n (v). For example, if the potential shape is 
f(x) — — sech 2 (a;), then the eigenvalues as functions of the coupling v are 
given 1 by the formula 



F n (v) 



l ~l 2 



n = 0,1,2,..., v>n(n + l). 

(1.1) 

These energy graphs are illustrated in Fig.(l). 

The problem we now consider here is the inverse of this, namely F — > f. 
We call this problem 'geometric spectral inversion'. It must at once be 
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E = F(v) for V(x) = -v sech A 2(x) 

1 | — 
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v 

Fig. 1. Discrete eigenvalues of H = — A — v sech 2 (x) as functions of the coupling v > 0. 

distinguished both from inverse scattering theory 2-6 and, more specifically, 
from the 'inverse problem in the coupling constant' discussed, for example, 
by Chadan et al. 2,7 ~ 10 In this latter problem, the discrete part of the 'input 
data' is a set {vi} of values of the coupling constant that all yield the 
identical energy eigenvalue E. The index i might typically represent the 
number of nodes in the corresponding eigenfunction. In contrast, for the 
problem discussed in the present paper, i is kept fixed and the input data 
is the graph (F(v), v), where the coupling parameter has any value v > v c , 
and v c is the critical value of v for the support of a discrete eigenvalue with i 
nodes. Geometric spectral inversion has been discussed in a series of earlier 
papers. n ~ 15 ' 24 Here we shall mainly discuss the bottom of the spectrum 
i = 0. However, on the basis of results we have obtained for the inversion 
IWKB of the WKB approximation, 13 there is good reason to expect that 
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inversion is possible starting from any discrete eigenvalue trajectory F n (v), 
n > 0. In fact, perhaps not surprisingly, IWKB yields better results starting 
from higher trajectories; moreover, they become asymptotically exact as the 
eigenvalue index is increased without limit. 

After recalling some basic results in the remainder of this introduction, 
in sections 2-4 we discuss a constructive algorithm for generating f(x) from 
F(v). In section-5 we apply this method to some examples. In sections 6 
and 7, we briefly describe an established geometric spectral theory 18-23 that 
enables us in section 8 to effect 'functional inversion'. 24 This latter method 
allows us to start with a seed potential and from this reconstruct 

f(x) by means of a sequence of functional operations. This latter method 
is applied to some examples in sections 9 and 10. 

By making suitable assumptions concerning the class of potential 
shapes, general theoretical progress has already been made with this inver- 
sion problem. 11,12 The assumptions that we retain for most of the present 
paper are that f(x) is symmetric, monotone increasing for x > 0, and 
bounded below: consequently the minimum value is /(0). We assume that 
our spectral data, the energy trajectory F(v), derives from a potential shape 
f(x) with these features. We have discussed 16 how two potential shapes /i 
and fi can cross over and still preserve spectral ordering F\ < F2. It is 
known 12 that lowest point /(0) of / is given by the limit 

/(0) = lim ^M. (1.2) 

v— >oo v 

We have proved 11 that a potential shape / has a finite flat portion (f'(x) — 
0) in its graph starting at x = if and only if the mean kinetic energy is 
bounded. That is to say, s — F(v) — vF'(v) < K, for some positive number 
K. More specifically, the size b of this patch can be estimated from F by 
means of the inequality: 

s<K => f(x) = f(0), \x\<b, and b=^K~i. (1.3) 

The monotonicity of the potential, which allows us to prove results like this, 
also yields the 
Concentration Lemma 11 

q(v) = J ip 2 {x,v)dx > f j^~_ F f ^ -> 1, as^tx), (1.4) 

where tp(x, v) is the normalized eigenfunction satisfying Hip = F{v)ip. More 
importantly, perhaps, if F(v) derives from a symmetric monotone poten- 
tial shape / which is bounded below, then / is uniquely determined. 12 
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The significance of this result can be appreciated more clearly upon con- 
sideration of an example. Suppose the bottom of the spectrum of H is 
given by F(v) = y/v, what is f{x)l It is well known, of course, that 
f(x) = x 2 — > F (v) = y/v: but are there any other potential pre- images 
of this spectral function F a (v)7 Are scaling arguments reversible? A pos- 
sible source of initial disquiet for anyone who ponders such questions is 
the uncountable number of (unsymmetric) perturbations 17 of the harmonic 
oscillator f(x) — x 2 there are, all of which have the identical spectrum 
E n = (2n + 1), n = 0, 1, 2, . . . , to that of the unperturbed oscillator (with 
coupling v = 1). 

If, in addition to symmetry and monotonicity, we also assume that a 
potential shape ,f a {x) vanishes at infinity and that f a (x) has area, then a 
given trajectory function F a (v) corresponding to f a (x) can be 'scaled' 12 to 
a standard form in which the new function F(v) — aF a ((3v) corresponds 
to a potential shape f(x) with area —2 and minimum value /(0) = — 1. 
Thus square- well potentials, which of course are completely determined 
by depth and area, are immediately invertible; moreover it is known that, 
amongst all standard potentials, the square-well is 'extremal' for it has the 
lowest possible energy trajectory. In Ref. 12 an approximate variational in- 
version method is developed; it is also demonstrated constructively that all 
separable potentials are invertible. However, these results and additional 
constraints are not used in the present paper. When a potential has area 
2A, we first assumed, during our early attempts at numerical inversion, 
that it would be very useful to determine A from F(v) and then appropri- 
ately constrain the inversion process. However, the area constraint did not 
turn out to be helpful. Thus the numerical method we have established for 
constructing f(x) from F(v) does not depend on use of this constraint, and 
is therefore not limited to the reconstruction of potentials which vanish at 
infinity and have area. 

2. Constructive inversion 

Much of numerical analysis assumes that errors arising from arithmetic 
computations or from the computation of elementary functions is negligibly 
small. The errors usually studied in depth are those that arise from the 
discrete representation of continuous objects such as functions, or from 
operations on them, such as derivatives or integrals. In this paper we shall 
take this separation of numerical problems to a higher level. We shall assume 
that we have a numerical method for solving the eigenvalue problem in 
the forward direction f(x) — > F(v) that is reliable and may be considered 
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for our purposes to be essentially error free. Our main emphasis will be 
on the design of an effective algorithm for the inverse problem assuming 
that the forward problem is numerically soluble. The forward problem is 
essential to our methods because we shall need to know not only the given 
exact energy trajectory F(v) but also, at each stage of the reconstruction, 
what eigenvalue a partly reconstructed potential generates. This line of 
thought immediately indicates that we shall also need a way of temporarily 
extrapolating a partly reconstructed potential to all x. 

Our constructive inversion algorithm hinges on the assumed symmetry 
and monotonicity of f(x). This allows us to start the reconstruction of 
f(x) at x = 0, and sequentially increase x. In Section (2) it is shown how 
numerical estimates can be made for the shape of the potential near x = 0, 
that is for x < 6, where 6 is a parameter of the algorithm. In Section (3) 
we explore the implications of the potential's monotonicity for the 'tail' of 
the wave function. In Section (4) we establish a numerical representation 
for the form of the unknown potential for x > b and construct our inversion 
algorithm. In Section (5) the algorithm is applied to three test problems. 

2.1. The reconstruction of f(x) near x = 0. 

Since the energy trajectory F(v) which we are given is assumed to arise 
from a symmetric monotone potential, and since the spectrum generated 
by the potential is invariant under shifts along the x-axis, we may assume 
without loss of generality that the minimum value of the potential occurs 
at x = 0. We now investigate the behaviour of F(v) 7 either analytically or 
numerically, for large values of v. The purpose is to establish a value for 
the starting point x = b > of our inversion algorithm and the shape of 
the potential in the interval x € [0, b}. First of all, the minimum value /(0) 
of the potential is provided by the limit (1.2). Now, if the mean kinetic 
energy s — (ip, — A^>) = F(v) — vF'(v) is found to be bounded above by a 
positive number K, then we know 11 that the potential shape f(x) satisfies 
f(x) = /(0), x e [0,6], where b is given by (1.3). In this case we have a 
value for 6 and also the shape f(x) inside the interval [0, 6]. 

If the mean potential energy s is (or appears numerically to be) un- 
bounded, then we adopt another strategy: we model f(x) as a shifted power 
potential near x — 0. Since we never know f(x) exactly, we shall need an- 
other symbol for the approximation we are currently using for f(x). We 
choose this to be g(x) and we suppose that the bottom of the spectrum of 
—A + vg(x) is given by G(v). The goal is to adjust g(x) until G(v) is close 
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to the given F(v). Thus we write 

f(x)ng(x) = f(0)+Ax q , xG[0,b]. (2.1) 

Therefore we have three positive parameters to determine, b, A, and q. We 
first suppose that g(x) has the form (2.1) for all x > 0. We now choose a 
'large' value v\ of v. This is related to the later choice of b by a bootstrap 
argument: the idea is that we choose v\ so large that the turning point 
determined by 

il>xx{x, vJMx, wi) = Vl f(x) - F(v) = (2.2) 

is equal to b. The concentration lemma guarantees that this is possible. By 
scaling arguments we have 

G(v) = f(0)v + E(q)(vA)&, (2.3) 

where E{q) is the bottom of the spectrum of the pure-power Hamiltonian 
-A + \x\ q . We now 'fit' G{v) to F(v) by the equations G(vi) = F(vi) and 
G{2v\) = F(2v\) which yield the estimate for q given by 

2 log(F(2^) - 2^/(0)) - log(F(^) - m/(0)) 
V= 2-T-q = ioi(2) • (2 ' 4) 

Thus A is given by 

A = ((F( Vl ) - «i/(0))/£7(g))^ (2.5) 

We choose 6 to be equal to the turning point corresponding to the model 
potential g(x) with the smaller value of v, that is to say so that f(0)+Ab q = 
F(vi)/vi, or 

b= ( F w-*m)\ ( 2 .5) 



Av\ 

Thus we have determined the three parameters which define the potential 
model g{x) for x G [—b, b]. 

3. The tail of the wavefunction 

Let us suppose that the ground-state wave function is ip(x,v). Thus the 
turning point tp xx (x, v) = occurs for a given v when 

x = x t (v) = r\R{v)), R(v)=(^y (3.1) 

The concentration lemma (1.4) quantifies the tendency of the wave function 
to become, as the coupling v is increased, progressively more concentrated 
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on the patch [— c, c], where x = c is the point (perhaps zero) where f(x) 
first starts to increase. This allows us to think in terms of the wave function 
having a 'tail'. We think of a symmetric potential as having been determined 
from x = up to the current point x. The question we now ask is: what 
value of v should wc use to determine how f(x), or, more particularly, 
our approximation g(x) for f(x), continues beyond the current point. We 
have found that a good choice is to choose v so that the turning point 
Xt{v) = x/2, or some other similar fixed fraction a < 1 of the current x 
value. The algorithm seems to be insensitive to this choice. Since g(x) has 
been constructed up to the current point, and F(v) is known, the value 
of v required follows by inverting (3.1). It has been proved 12 that R(v) 
is monotone and therefore invertible. Hence we have the following general 
recipe for v : 

v = R-\g(ax)), a = \. (3.2) 

Since we can only determine Schrodinger eigenvalues of H = — A + vg{x) if 
the potential is defined for all x, we must have a policy about temporarily 
extending g(x). We have tried many possibilities and found the simplest 
and most effective method is to extend g(x) in a straight line, with slope 
to be determined. 

4. Constructive inversion algorithm 

We must first define the 'current' approximation g(x) for the potential f(x) 
sought. For values of x less than b, g(x) is defined either as the horizontal 
line f(x) = /(0) or as the shifted power potential (2.1). For values of 
x greater than b, the x-axis is divided into steps of length h. Thus the 
'current' value of x would be of the form x = x k = b + kh, where k is a 
positive integer. The idea is that g(xk) is determined sequentially and g{x) 
is interpolated linearly between the Xk points. We suppose that {g{xk)} 
have already been determined up to k and we need to find y = g(xk+i)- 
For x > Xk we let 

g(x) = g{x k ) + (y - g(x k )) — — — . (4.1) 

If, from a study of F(v), the underlying potential f(x) has been shown 12 to 
be bounded above, it is convenient to rescale F(v) so that it corresponds to a 
potential shape f(x) which vanishes at infinity. In this case it is slightly more 
efficient to modify (4.1) so that for large x the straight-line extrapolation 
of g(x) is 'cut' to zero instead of becoming positive. In either case we now 
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Fig. 2. Wc illustrate the ideas discussed in the text for the case of the sech-squarcd 
potential. The inset graph shows the sech-squared potential perturbed from x = x a by 
five straight line extensions; meanwhile the main graph shows the corresponding set of 
five wave functions which agree for < x < x a and then continue with different 'tails' 
dictated by the corresponding potential extensions. The value of the coupling v is the 
value that makes the turning point of the wave function occur at x = x a /2. This figure 
illustrates the sort of graphical study that has lead to the algorithm described in this 
paper. 

have for the current point Xk an approximate potential g(x) parameterized 
by the 'next' value y = g(xk+i)- The task of the inversion algorithm is 
simply to choose this value of y. 

Let us suppose that, for given values of k and y, the bottom of the 
spectrum of H = —A + vg(x) is given by G(v,k,y), then the inversion 
algorithm may be stated in the following succinct form in which a < 1 is a 
fixed parameter. Find y such that 



vg(ax k ) = F(v) = G(v, k, y); then g(x k+1 ) = y. (4.2) 
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The value of v is first chosen so that the turning point of the wave function 
generated by g occurs at axk', after this, the value of y is chosen so that 
G 'fits' F for this value of v. The value of the parameter a chosen for the 
examples discussed in section (5) below is a = |. The idea behind this 
choice can best be understood from a study of Figure (1): the value of the 
coupling v must be such that the current value of x for which y is sought is 
in the 'tail' of the corresponding wave function; that is to say, the turning 
point ax should be before x, but not too far away. Fortunately the inversion 
algorithm seems to be insensitive to the choice of a. 

5. Some examples 

The first example we consider is the unbounded potential whose shape f(x) 
and corresponding exact energy trajectory F(v) are given by the {/, F} pair 

f(x) = -l + \x\* i — > F(v) = -v + E{3/2)v$, (5.1) 

where £7(3/2) is the bottom of the spectrum of H = — A+ \x\i and has the 
approximate value E(3/2) w 1.001184. Applying the inversion algorithm to 
F(v) we obtain the reconstructed potential shown in Figure (3). We first 
set v\ = 10 4 and find that the initial shape is determined (as described in 
Section (2)) to be —1 + x 1 ' 5 for x < b = 0.072. For larger values of x the 
step size is chosen to be h = 0.05 and 40 iterations are performed by the 
inversion algorithm. The results are plotted as hexagons on top of the exact 
potential shape shown as a smooth curve. 

The following two examples are bounded potentials both having large-x 
limit zero, lowest point /(0) = — 1, and 'area' —2. The exponential poten- 
tial 1 - 26 has the {f,F} pair 

f(x) = -e-l x l «— > J' L (2vi)=0 = E = F{v) 1 (5.2) 

2|£7j 2 

where J' v (x) is the derivative of the Bessel function of the first kind of order 
v. 

For the sech-squared potential 1 we have 

2 

(5.3) 

In Figure (4) the two energy trajectories arc plotted. Since the two potentials 
have lowest value —1 and 'area' —2 it follows 12 that the corresponding 
trajectories both have the form F(v) s=s — v 2 for small v and they both satisfy 
the large-w limit lim^oo (F(v)/v) = —1. Thus the differences between the 



f(x) = — sech 2 (x) 



F (v) 
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f(x) 




X 



Fig. 3. Constructive inversion of the energy trajectory F(v) for the shifted power po- 

3 

tential f(x) = —1 + \x\? . For x < b = 0.072, the algorithm correctly generates the 
model f(x); for larger values of x, in steps of size h = 0.05, the hexagons indicate 
the reconstructed values for the potential f(x), shown exactly as a smooth curve. The 
unnormalized wave functions are also shown. 

potential shapes is somehow encoded in the fine differences between these 
two similar energy curves for intermediate values of v : it is the task of 
our inversion theory to decode this information and reveal the underlying 
potential shape. If we apply the inversion algorithm to these two problems 
we obtain the results shown in Figures (5) and (6). The parameters used 
are exactly the same as for the first problem described above. The time 
taken to perform the inversions was less than 20 seconds if we discount, in 
the case of the exponential potential, the extra time taken to compute F(v) 
itself. 

6. Kinetic potentials 

Geometry is involved with this problem because we deal with a family of 
operators depending on a continuous parameter v. This immediately leads 
to a family of spectral manifolds, and, more particularly, to the consider- 
ation of smooth transformations of potentials, and to the transformations 
which they in turn induce on the spectral manifolds. This is the environ- 
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Fig. 4. The ground-state energy trajectories F(v) for the exponential potential (E) and 
the sech-squarcd potential (S). For small v, F(v) —v ; for large v, lim^-joo (F(v)/v) = 
— 1. The shapes of the underlying potentials are buried in the details of F(v) for inter- 
mediate values of v. 

ment in which we are able to construct the following functional inversion 
sequence: 

f[n+l] = f Q fin]- 1 Q f[n] =J Q R [n] _ (g^ 

A kinetic potential is the constrained mean value of the potential shape 
f( s ) = (/); where the corresponding mean kinetic energy s = (—A) is 
held constant. It turns out that kinetic potentials may be obtained from 
the corresponding energy trajectory F by what is essentially a Legendre 
transformation 25 / <H> F given 20 by 

f(s) = F'(v), s = F(v)-vF'(v), 

and 

F(v)/v = f(s)-sf'(s), l/v = -f'{s). (6.2) 

These transformations are well defined because of the definite convexities 
of F and /; they complete the definition of the inversion sequence (1.2), up 
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Fig. 5. Constructive inversion of the energy trajectory F(v) for the exponential poten- 
tial f(x) = — exp(:r). For x < b = 0.048, the algorithm correctly generates the model 
f(x) = —1 + |x|; for larger values of x, in steps of size h = 0.05, the hexagons indicate 
the reconstructed values for the potential f(x), shown exactly as a smooth curve. The 
unnormalized wave functions are also shown. 

to the choice of a starting seed potential f^(x). They differ from Legendre 
transformations only because of our choice of signs. The choice has been 
made so that the eigenvalue can be written (exactly) in the semi-classical 
forms 

E = F{v) = mm{s + vf(s)} = min Ik^Mx) + vf(x)\ (6.3) 

s>0 x>0 I J 

where the kinetic- and potential-energy terms have the 'usual' signs. 
7. Envelope theory 

The term 'kinetic potential' is short for 'minimum mean iso-kinetic po- 
tential'. If the Hamiltonian is H = — A + vf(x), where f(x) is potential 
shape, and T>(H) C L 2 ($i) is the domain of H, then the ground-state kinetic 
potential f(s) — fo{s) is defined 19,20 by the expression 
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o 



f(x) 



Fig. 6. Constructive inversion of the energy trajectory F(v) for the sech-squared po- 
tential f(x) = — sech 2 (x). For x < b = 0.1, the algorithm correctly generates the model 
f(x) = — 1 + x 2 ; for larger values of x, in steps of size h = 0.05, the hexagons indicate 
the reconstructed values for the potential f(x), shown exactly as a smooth curve. The 
unnormalized wave functions are also shown. 



/»= Jn^ (7.1) 

(V,V)=i 

(i/>,-Ai/>)=s 

The extension of this definition to the higher discrete eigenvalues (for 
v sufficiently large) is straightforward 20 but not explicitly needed in the 
present paper. The idea is that the min-max computation of the discrete 
eigenvalues is carried out in two stages: in the first stage (7.1) the mean 
potential shape is found for each fixed value of the mean kinetic energy s; 
in the second and final stage we minimize over s. Thus we have arrive at 
the semi-classical expression which is the first equality of Eq.(1.4). It is well 
known that F(v) is concave (F"(v) < 0) and it follows immediately that 
f(s) is convex. More particularly, we have 11 



F"{v)j"{s) = ~. (7.2) 
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Thus, although kinetic potentials are defined by (7.1), the transforma- 
tions (6.2) may be used in practice to go back and forth between F and 

Kinetic potentials have been used to study smooth transformations of 
potentials and also linear combinations. The present work is an application 
of the first kind. Our goal is to devise a method of searching for a trans- 
formation g, which would convert the initial seed potential /'°'(a;) into the 
(unknown) goal f{x) — g(f^)- We shall summarize briefly how one pro- 
ceeds in the forward direction, to approximate F, if we know f(x). The K 
functions are then introduced, by a change of variable, so that the potential 
f(x) is exposed and can be extracted in a sequential inversion process. 

In the forward direction we assume that the lowest eigenvalue F^(v) 
of £f[°l = — A + w/[°](x) is known for all v > and we assume that f(x) is 
given; hence, since the potentials are symmetric and monotone for x > 0, 
we have defined the transformation function g. 'Tangential potentials' to 
g(/[°l) have the form a+bf^(x), where the coefficients a(t) and b(t) depend 
on the point of contact x = t of the tangential potential to the graph of 
f(x). Each one of these tangential potentials generates an energy trajectory 
of the form T{v) = av + F^(bv), and the envelope of this family (with 
respect to t) forms an approximation F A (v) to F(v). If the transformation 
g has definite convexity, then F A (v) will be either an upper or lower bound 
to F(v). It turns out 20 that all the calculations implied by this envelope 
approximation can be summarized nicely by kinetic potentials. Thus the 
whole procedure just described corresponds exactly to the expression: 

f*f A = g°f m , (7.3) 

with w being replaced by an inequality in case g has definite convexity. Once 
we have an approximation f A , we immediately recover the corresponding 
energy trajectory F A from the general minimization formula (6.3). 

The formulation that reveals the potential shape is obtained when we 
use x instead of s as the minimization parameter. We achieve this by the 
following general definition of x and of the K function associated with / : 

f(x) = f(s), KW{x) = r\f(x)). (7.4) 

The monotonicity of f(x) and of / guarantee that x and K are well defined. 
Since / _1 (/) is a convex function of /, the second equality in (1.4) imme- 
diately follows. 19 In terms of K the envelope approximation (2.3) becomes 
simply 

KWnKi'™]. (7.5) 
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Thus the envelope approximation involves the use of an approximate K 
function that no longer depends on /, and there is now the possibility that 
we can invert (1.4) to extract an approximation for the potential shape. 

We end this summary by listing some specific results that we shall need. 
First of all, the kinetic potentials and K functions obey 18,19 the following 
elementary shift and scaling laws: 

f(x) -> a + bf(x/t) => {/» -+ a + bf(st 2 ), KW{x) -+ ^K^ (|) 

(7.6) 

Pure power potentials are important examples which have the following 
formulas: 

f(x)=\x\«^i^f(s)= gy, K(x)=^yy n.i) 

where, if the bottom of the spectrum of —A + \x\ q is E(q), then the P 
numbers are given 19 by the following expressions with n = : 

i i 

Pn(q) = \E n {q)\^ 

We have allowed for q < and for higher eigenvalues since the formulas 
are essentially the same. The P n (q) as functions of q are interesting in 
themselves: 19 they have been proved to be monotone increasing, they are 
probably concave, and P n (0) corresponds exactly to the log potential. By 
contrast the E n (q) are not so smooth: for example, they have infinite slopes 
at q = 0. But this is another story. An important observation is that the 
K functions for the pure powers are all of the form (P(q)/x) 2 and they 
are invariant with respect to both potential shifts and multipliers: thus 
a + b\x\ 9 has the same K function as does \x\ q . For the harmonic oscillator 

P„(2) = (n + i) 2 , n = 0, 1, 2, Other specific examples may be found 

in the references cited. 

The last formulas we shall need are those for the ground state of the 
sech-squared potential: 

f{x) - -sech 2 (x) => [fis) = — ^ , Kix) = sinh- 2 ^)) . 

(7.9) 

8. Functional inversion 

The inversion sequence (6.1) is based on the following idea. The goal is to 
find a transformation g so that / = g o We choose a seed f^°\ but, 



2 




M 


2 + q 




2 + q 



q^O. 



(7.8) 
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of course, / is unknown. In so far as the envelope approximation with 

as a basis is 'good', then an approximation </M for g would be given by 

f = g m o Thus we have 

gn g m=fofM-\ (8.1) 
Applying this approximate transformation to the seed we find: 

/ w /W = gW o /P] = / o /I"]" 1 o /P] = / o Jft"]. (8.2) 

We now use /W as the basis for another envelope approximation, and, 
by repetition, we have the ansatz (1.2), that is to say 

f[n+l] = j Q fin]- 1 Q f[n] = j Q R [n] ^ 

A useful practical device is to invert the second expression for F given 
in (1.4) to obtain 

K^(x)=max{F(v)-vf(x)\. (8.4) 

v>0 

The concavity of F(v) explains the max in this inversion, which, as it 
stands, is exact. In a situation where / is unknown, we have / on both sides 
and nothing can be done with this formal result. However, in the inversion 
sequence which we are considering, (3.4) is extremely useful. If we re-write 
(3.4) for stage [n] of the inversion sequence it becomes: 

(x) = max (fW (v) - w/M (a:)) . (8.5) 

v>0 I J 

In this application, the current potential shape and consequently 
can be found (by shooting methods) for each value of v. The mini- 
mization can then be performed even without differentiation (for example, 
by using a Fibonacci search) and this is a much more effective method for 
K [n] = jln]- Q fin] than filing /W( s ), finding the functional inverse, and 
applying the result to /M. 

9. Functional inversion for pure powers 

We now treat the case of pure-power potentials given by 
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f(x) = A + B\x\\ q>0, (9.1) 

where A and B > are arbitrary and fixed. We shall prove that, starting 
from another pure power as a seed, the inversion sequence converges in 
just two steps. The exact energy trajectory F(v) for the potential (9.1) is 
assumed known. Hence, so is the exact kinetic potential given by (7.7) and 
the general scaling rule (7.6), that is to say 

f{s) = A + B(^y . (9.2) 
We now suppose that a pure power is also used as a seed, thus we have 



fM(x)=a + b\x\ p => KM(x)= , (9.3) 

where the parameters a, b > 0, p > are arbitrary and fixed. The 
first step of the inversion (6.3) therefore yields 

= (/o*M) {x) =A + B (^My. (9.4) 

The approximate potential f^(x) now has the correct x power depen- 
dence but has the wrong multiplying factor. Because of the invariance of 
the K functions to multipliers, this error is completely corrected at the next 
step, yielding: 



K[1]{x ^ = \x) ^ f m {x)={foK^{x)=A + B\x\ q . (9.5) 

This establishes our claim that power potentials are inverted without 
error in exactly two steps. 

The implications of this result are a little wider than one might first 
suspect. If the potential that is being reconstructed has the asymptotic 
form of a pure power for small or large x, say, then we know that the 
inversion sequence will very quickly produce an accurate approximation for 
that part of the potential shape. More generally, since the first step of the 
inversion process involves the construction of K^°\ the general invariance 
property K^ a+b ^ = given in (7.6) means that the seed potential 
may be chosen without special consideration to gross features of / already 
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arrived at by other methods. For example, the area (if the potential has 
area), or the starting value /(0) need not be incorporated in say, by 
adjusting a and b. 



10. A more general example 

We consider the problem of reconstructing the sech-squared potential 
f(x) = -secli (x). We assume that the corresponding exact energy trajec- 
tory F(v) and, consequently, the kinetic potential f(s) are known. Thus 19 
from the potential shape f(x) = — sech 2 (a;) we have: 



F(v) = 



1\ 3 1 

y+ 4 "2 



m = - ( %r- } (io.i) 

(s + s 2 j 2 + s 



. The seed is essentially x 2 , but we use a scaled version of this for the 
purpose of illustration in Fig. (7). Thus we have 

/ [0] = -l + g => K [0] (*) = ^ (10-2) 
This potential generates the exact eigenvalue 

FW ^ = - v+ {h) i (10 ' 3) 
, which, like the potential itself, is very different from that of the goal. After 
the first iteration we obtain 

/W(x) = / ( K M(x)) = - — -2- (10.4) 

\ ' l + (l + 4a; 2 )2 

A graph of this potential is shown as /l in Fig. (7). In order to continue 
analytically we would need to solve the problem with Hamiltonian i/W = 
—A + u/M(ar) exactly to find an expression for F^(v). We know no way 
of doing this. However, it can be done numerically, with the aid of the 
inversion formula (8.5) for K. The first 5 iterations shown in Fig. (7) suggest 
convergence of the functional sequence. 



11. Conclusion 

We have discussed two methods that may be used to reconstruct the po- 
tential f(r) in a Schrodinger operator H = — A + vf(r) if an eigenvalue 
curve E = F n (v) is known as a function of the coupling v. Such functions 
as F(v) are met, for example, as one-body approximations to certain iV- 
body systems, such as atoms. In addition to the inversion of physical data 



March 8, 2013 11:58 



WSPC - Proceedings Trim Size: 9in x 6in rcmpl3hall 



19 




Fig. 7. The energy trajectory F for the sech-squared potential f(x) = — sech 2 (:r) is 
approximately inverted starting from the seed f^(x) = —1 + x 2 /20. The first step can 
be completed analytically yielding fl = /M (x) = — 2/{l + Vl + 4a; 2 }. Four more steps 
{fk = f lk] }U 2 of the inversion sequence approaching / are performed numerically. 



to infer the structure of the underlying system, the situation also presents 
a rather fascinating mathematical problem. Work is under way extending 
these results and methods to a wider class of problems, including those with 
singular potentials. 
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